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Abstract 

A boundary control problem for the pure Gahn-Hilliard equations with possi¬ 
bly singular potentials and dynamic boundary conditions is studied and first-order 
necessary conditions for optimality are proved. 
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1 Introduction 

The simplest form of the Cahn-Hilliard equation (see IHElllla]) reads as follows 

dtU — Aw = 0 and w = —Ay + f'{y) in P x (0, T), (1.1) 
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PRIN Grant 2010A2TFX2 “Calculus of Variations” and the GNAMPA (Gruppo Nazionale per I’Analisi 
Matematica, la Probabilita e le loro Applicazioni) of INdAM (Istituto Nazionale di Alta Matematica). 
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where 12 is the domain where the evolution takes place, and y and w denote the order 
parameter and the chemical potential, respectively. Moreover, /' represents the derivative 
of a double well potential /, and typical and important examples are the following 

freg{r) = \{T^-l)\ TGR (1.2) 

/«og(p) = ((1 + p) ln(l + r) + (1 - r) ln(l - r)) - cr^ , r G (-1,1), (1.3) 

where c > 0 in the latter is large enough in order that fiog be nonconvex. The potentials 
fll.2p and fll.3p are usually called the classical regular potential and the logarithmic double¬ 
well potential, respectively. 

The present paper is devoted to the study of the control problem described below 
for the initial-boundary value problem obtained by complementing fll.ip with an initial 
condition like i/(0) = yo and the following boundary conditions 

dnW = 0 and dny + dtyr - ^ryr + frivr) = ur on T x (0, T) (1.4) 

where T is the boundary of 12. The former is very common in the literature and preserves 
mass conservation, i.e., it implies that the space integral of y is constant in time. The 
latter is an evolution equation for the trace yr of the order parameter on the boundary, 
and the normal derivative dny and Mr act as forcing terms. This condition enters the 
class of the so-called dynamic boundary conditions that have been widely used in the 
literature in the last twenty years, say: in particular, the study of dynamic boundary 
conditions with Cahn-Hilliard type equations has been taken up by some authors (let us 
quote [51 El d CSl ca 121] and also refer to the recent contribution [S] in which also a 
forced mass constraint on the boundary is considered). 

The dynamic boundary condition in fll.dp contains the Laplace-Beltrami operator Ap 
and a nonlinearity /f which is analogous to f but is now acting on the boundary values Mr- 
Even though some of our results hold under weaker hypotheses, we assume from the very 
beginning that /' and /f have the same domain 2). The main assumption we make 
is a compatibility condition between these nonlinearities. Namely, we suppose that /f 
dominates f in the following sense: 

\f'{r)\ <r]\f^{r)\+C (1.5) 

for some positive constants y and C and for every r G D. This condition, earlier introduced 
in [1] in relation with the Allen-Cahn equation with dynamic boundary conditions (see 
also d), is then used in [9| (as well as in [6| and cni) to deal with the Cahn-Hilliard 
system. This complements [T3|, where some kind of an opposite inequality is assumed. 

As just said, this paper deals with a control problem for the state system described 
above, the control being the source term Mr that appears in the dynamic boundary con¬ 
dition fll.4p . Namely, the problem we want to address consists in minimizing a proper 
cost functional depending on both the control Mr and the associate state {y,yr). Among 
several possibilities, we choose the cost functional 

3{y,yr,ur) := y \\y - zq \\ 12 (^ q ^ + y ||i/r - ^E||i 2 (s) + y WurWh^^s) > (1-6) 

where the functions zq , zj : and the nonnegative constants are given. The control 

problem then consists in minimizing fll.bp subject to the state system and to the constraint 
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Mr £ U^ad, where the control box Had is given by 

Uad:= {ureH\0,T-,Hr)nL^{E)-. 

Ur ,min < Mr < Rr,max a-e. on S, ll^tMr< Mq} (1.7) 

for some given functions Mr, min, Rr.max ^ and some prescribed positive constant Mq. 

Of course, the control box Had must be nonempty and this is guaranteed if, for instance, 
at least one of Mr,min or Mr,max is in 77^(0, T; Hr) and its time derivative satishes the above 
L^(S) bound. 

This paper is a follow-up of the recent contributions [9] and [10] already mentioned. 
They deal with a similar system and a similar control problem. The paper [9| contains a 
number of results on the state system obtained by considering 


m ^ rdifi - Ay + f'(y) 


( 1 . 8 ) 


in place of the second condition in fll.ip . In fll.Sp . r is a nonnegative parameter and 
the case r > 0 coupled with the hrst equation in fll.ll) yields the well-known viscous 
Cahn-Hilliard equation (in contrast, we term (11.11) the pure Cahn-Hilliard system). More 
precisely, existence, uniqueness and regularity results are proved in [9] for general poten¬ 
tials that include fll.2l) - fll.3l) . and are valid for both the viscous and pure cases, i.e., by 
assuming just r > 0. Moreover, if r > 0, further regularity and properties of the solution 
are ensured. These results are then used in ra. where the boundary control problem 
associated to a cost functional that generalizes fll.61) is addressed and both the existence 
of an optimal control and hrst-order necessary conditions for optimality are proved and 
expressed in terms of the solution of a proper adjoint problem. 

In fact, recently Cahn-Hilliard systems have been rather investigated from the view¬ 
point of optimal control. In this connection, we refer to [151 ESI ET] and to [25l 126] which 
deal with the convective Cahn-Hilliard equation; the case with a nonlocal potential is 
studied in izni. There also exist contributions addressing some discretized versions of 
general Cahn-Hilliard systems, cf. [TBl E2] • However, about the optimal control of vis¬ 
cous or non-viscous Cahn-Hilliard systems with dynamic boundary conditions of the form 
fll.41) . we only know of the papers [10] and [6] dealing with the viscous case; to the best 
of our knowledge, the present contribution is the hrst paper treating the optimal control 
of the pure Cahn-Hilliard system with dynamic boundary conditions. 

The technique used in our approach essentially consists in starting from the known 
results for r > 0 and then letting the parameter r tend to zero. In doing that, we 
use some of the ideas of [7] and [6], which deal with the Allen-Cahn and the viscous 
Cahn-Hilliard equations, respectively, and address similar control problems related to the 
nondifferentiable double obstacle potential by seeing it as a limit of logarithmic double¬ 
well potentials. 

The paper is organized as follows. In the next section, we list our assumptions, state 
the problem in a precise form and present our results. The corresponding proofs are given 
in the last section. 
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2 Statement of the problem and results 


In this section, we describe the problem nnder study and give an outline of our results. 
As in the Introduction, is the body where the evolution takes place. We assume C 
to be open, bounded, connected, and smooth, and we write |f2| for its Lebesgue measure. 
Moreover, F, Vr and Ar stand for the boundary of 14, the outward normal derivative, 
the surface gradient and the Laplace-Beltrami operator, respectively. Given a hnite hnal 
time T > 0, we set for convenience 

Qt := 14 X (0,4) and := F x (0,4) for every 4 G (0,T] (2.1) 

Q '.= Qj ', and B := Bj’. (2-2) 


Now, we specify the assumptions on the structure of our system. Even though some of 
the results we quote hold under rather mild hypotheses, we give a list of assumptions that 
implies the whole set of conditions required in [9]. We assume that 


—oo < r_ < 0 < r+ < +00 
/, /r : (r_,r+) —)■ [0, +cxo) are functions 
/(O) = /r(0) = 0 and /" and /f' are bounded from below 
|/'(r)| < 7 ] |/f(r)| + C for some ? 7 , G > 0 and every r G (r_,r+) 

and 


lim f'{r) = lim /f(r) = —oo 

rN^r— F\,r_ 


lim f'{r) = lim /f (r) = +cxo 


(2.3) 

(2.4) 

(2.5) 

(2.6) 

(2.7) 


We note that fl2.3l) - fl2.7l) imply the possibility of splitting /' as /' = /3 + tt, where /3 is a 
monotone function that diverges at r± and vr is a perturbation with a bounded derivative. 
Moreover, the same is true for /r, so that the assumptions of [9] are satished. Furthermore, 
the choices / = freg and / = fiog corresponding to fll.2p and fll.31) are allowed. 

Next, in order to simplify notations, we set 

V:=H\n), H:=L‘^{n), Hr := L‘^{T) and Vr := H^T) (2.8) 

V := {(n,nr) G B x W : Pr = P|p} and !K := H x Hr (2.9) 

and endow these spaces with their natural norms. If X is any Banach space, then 

II • \\x and X* denote its norm and its dual space, respectively. Furthermore, the symbol 

(•, •) usually stands for the duality pairing between V* and V itself and the similar no¬ 
tation (•, • )r refers to the spaces and Vr. In the following, it is understood that H 
is identihed with H* and thus embedded in V* in the usual way, i.e., such that we have 
{u,v) = {u,v) with the inner product {■,■) of H, for every u E H and v E V. Thus, we 
introduce the Hilbert triplet (V, H, V*) and analogously behave with the boundary spaces 
Vr, Hr and V^. Finally, ifuEV* and u E L^(0,T; V*), we dehne their generalized mean 
values G M and yP E L^{0, T) by setting 

:= {u,l) and m*^(4) := (w(4))^ for a.a. 4 G (0,T). (2-10) 


Clearly, the relations in fl2.10p give the usual mean values when applied to elements of H 
or H{f),T-,H). 
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At this point, we can describe the state problem. For the data, we assume that 


and e H\T) (2.11) 

r_ < |/o(x) < r+ for every x G (2-12) 

ur&H\fd,T-Hr). (2.13) 

We look for a triplet (y, i/r, w) satisfying 

y G T; V*) C L°°(0, T; V) C T; H^{n)) (2.14) 

yr G H\0, T; Hr) H L-(0, T; W) C ^^(0, T; H^iT)) (2.15) 

yr{t) = y{t)\^ for a.a. t G (0,T) (2.16) 

weL^{0,T;V), (2.17) 

as well as, for almost every t G (0,T), the variational equations 


{dty{t) v)+ ■ Vn = 0 for every n G C 

Jn 

I w{t)v= I Vy{t)-Vv+ I dtyr{t)vr+ [ 'Vryr{t) ■ VrVr 
Jn Jn Jr Jr 

+ [ fivit)) 'V+ [ (friyrit)) - ur{t)) vr for every (n, wr) e V 
Jn Jr 


»(•)) = yo ■ 


(2.18) 


(2.19) 

( 2 . 20 ) 


Thus, we require that the state variables satisfy the variational counterpart of the problem 
described in the Introduction in a strong form. We note that an equivalent formulation 
of fl2T8ll - flCTll is given by 


{dty{t) v{t)) dt + / Vtc • Vn = 0 

Jq 

wv= Vy • Vn + / dtyrvr+ / Vr|/r • VrUr 


( 2 . 21 ) 


'Q 


>Q 


+ / f{y) v+ {friyr) - Ur) vr 


( 2 . 22 ) 


’Q 


for every v G L^(0,T;C) and every (u,Pr) G L^(0,T;V), respectively. 

Besides, we consider the analogous state system with viscosity. Namely, for r > 0 we 
replace fl2.19p by 

/ w{t)v = T / dty{t)v+ / Vy{t)-Vv+ / dtyr{t)vr+ / Vr|/r(t) ■ Vrur 
Jn Jn Jn Jr Jr 

+ [ f'{y{t))v+ [ {fr{yr{t)) - Ur{t)) Vr for every (u, up) e V (2.23) 

Jn Jr 

in the above system. We notice that a variational equation equivalent to fl2.23p is given 
by the analogue of (I2.22p . i.e.. 


wv = T dtyv+ / Vy-Vv+ / dtyrvr+ / Vpi/r ■ Vpfr 


'Q 


'Q 


'Q 


+ / f{y)v+ / {friyr) - Ur) Vr for every (u, up) e L2(0, T; V). (2.24) 


’Q 
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As far as existence, uniqueness, regularity and continuous dependence are concerned, 
we directly refer to [5]. From 0 Thms. 2.2 and 2.3] (where V has a slightly different 
meaning with respect to the present paper), we have the following results: 

Theorem 2 . 1 . Assume fl2.3p - fl2.7p and fl2.lip - fl2.13p . Then, there exists a unique triplet 
{y,yr,w) satisfying fl2.14p - fl2.17p and solving fl2.18p - fl2.20p . 

Theorem 2 . 2 . Assume fl2.3p - fl2.7p and fl2.lip - fl2.13p . Then, for every r > 0, there exists 
a unique triplet {y'^,yf,w'^) satisfying fl2.14p - fl2.17p and solving fl2.18p . fl2.20p and fl2.23p . 
Moreover, this solution satisfies dty'^ G L^(0,T;i7) and the estimate 

||2/’'||Rl(0,r;U*)nL°°(0,T;y)nL2(0,r;R2(Q)) 

+ 11 l/f 11 (0,T;Rr)nL°° {0,T-,Vr)nL^ (F)) 

+ II'W^’^I|l2(0,T;U) + \\f'{y^)\\L^{0,T-,H) + \\ friyr)\\ L^{0,T-,Hr) 

+ T^^‘^\\dty^\\L^{0,T;H) < C'o (2.25) 

holds true for some constant Co > 0 that depends only on fl, T, the shape of the non- 
linearities f and fr, and the norms ||(|/o,|/o|r)||v^ ||/'(2 /o)||l1(o); ||/f(l/o|r) 
||wr||L2(0,T;Rr)- 

In fact, if the data are more regular, in particular, if Mr ^ -^^(0) T ; Hr) nL°°(E), then 
the solution [y'^,yf,w'^) is even smoother (see 0 Thms. 2.4 and 2.6]) and, specihcally, it 
satishes 

ph < a.e. in Q (2.26) 

for some constants rh,r^ G (r_,r+) that depend on r, in addition. It follows that the 
functions f"{y'^) and fr{yf) (which will appear as coefficients in a linear system later on) 
are bounded. We also notice that the stability estimate fl2.25p is not explicitely written 
in 0. However, as the proof of the regularity fl2.14p - fl2.17p of the solution performed 
there relies on a priori estimates and compactness arguments and the dependence on r 
in the whole calculation of 0 is always made explicit, fl2.25p holds as well, and we stress 
that the corresponding constant Co does not depend on r. 

Once well-posedness for problem fl2.18l) - fl2.20l) is established, we can address the cor¬ 
responding control problem. As in the Introduction, given two functions 

zq G L2(Q) and Zs e L‘^{T) (2.27) 

and three nonnegative constants bg, 6s, 6o, we set 

3{y,yr,ur) ■■= lb - ^qIIl 2 (q) + Y Ibr - ^s||i 2 (s) + ^ Ibr|li 2 (s) (2.28) 

for, say, y G L^(0,T;i7), yr G L^(0,T;i7r) and Mr £ -^^(^), and consider the problem 
of minimizing the cost functional fl2.28p subject to the constraint Mp G Had, where the 
control box Had is given by 

Had:= {Mr G 77^(0, T; i7r) n L°°(E) : 

Rr,min ^ Rr ^ Rr,max a.e. on B, ||c?tMr||L2(E) ^ Afoj" (2.29) 
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and to the state system fl2.18l) - fl2.20l) . We simply assume that 

Mq > 0, Mr,min, Mr,max ^ L°°(E) and Had is nonempty. (2.30) 

Besides, we consider the analogous control problem of minimizing the cost functional 
fl2.28p subject to the constraint Mr G hlad and to the state system 02.181) . O2.20p and 02.231) . 
From dQi Thm. 2.3], we have the following result. 

Theorem 2.3. Assume 02.3p - 02.7p and 02.1ip - 02.13p . and let 3 and Had be defined by 
02.28P and 02.29P under the assumptions 02.27P and O2.30p . Then, for every r > 0, there 
exists Mf G had such that 

3{y^, Vr, Mf) < 3{y'", yf, Mr) for every Mr G had , (2.31) 

where y'^, yfi, y'^ and yf are the components of the solutions ,yf ,w'^) and {y'^,yf,w'^) 
to the state system 02.18p . O2.20p and 02.231) corresponding to the controls uf and Mr, 
respectively. 

In [TU] hrst-order necessary conditions are obtained in terms of the solution to a proper 
adjoint system. More precisely, just the case r = 1 is considered there. However, by going 
through the paper with some care, one easily reconstructs the version of the adjoint system 
corresponding to an arbitrary r > 0. Even though the adjoint problem considered in Hi 
involves a triplet (p"^, qf) as an adjoint state, only the third component qf enters the 
necessary condition for optimality. On the other hand, q'^ and qf are strictly related to 
each other. Hence, we mention the result that deals with the pair (g"^, gf). To this end, we 
recall a tool, the generalized Neumann problem solver N, that is often used in connection 
with the Cahn-Hilliard equations. With the notation for the mean value introduced 
in fl2.10p . we dehne 

domTsf := {m* G : m)]^ = 0} and N : domN {v E V : m^ = 0} (2.32) 

by setting, for m* G domN, 

Nm* G V, (Nm*)*^ = 0, and / VNm* • Vz = (m*, 2 ;) for every z eV. (2.33) 

Jn 

Thus, Nm* is the solution v to the generalized Neumann problem for —A with datum m* 
that satishes = 0. Indeed, if m* G H, the above variational equation means that 
—ANm* = M* and = 0. As 12 is bounded, smooth, and connected, it turns out that 

fl2.33p yields a well-dehned isomorphism. Moreover, we have 

(m*,Nm*) = (m*,Nm*) = / (VNm*) • (VNm*) for m*,m* G domN, (2.34) 

Jn 

whence also 

2{dtvfit),JJvfit)) = — J |VNM*(t)r = ^ ||M*(f)||* for a.a. f G (0,T) (2.35) 

for every m* G i/^(0,T; V*) satisfying (m*)^ = 0 a.e. in (0,T), where we have set 


:= ||VNM*||i^ for M* G V*. 


(2.36) 
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One easily sees that || • ||* is a norm in V* which is equivalent to the usual dual norm. 
Furthermore, we introduce the spaces and Vq by setting 

^Kn:={(p,Pr) eF{: u^ = 0} and Vn:=F£f,nV, (2.37) 

and endow them with their natural topologies as subspaces of dC and V, respectively. As 
in m Thms. 2.5 and 5.4], we have the following result. 

Theorem 2.4. Assume 

A G Ar e ifQ G L^{Q) and (ps e L\E). (2.38) 

Then, for every t > 0, there exists a unique pair (g"^, qf) satisfying the regularity conditions 
q^ e and qf & H\Q,T-Hr) A H‘^{V)) (2.39) 

and solving the following adjoint problem: 


e Vf^ for every t G [0,T] 

- [ dt{'N{q^{t)) + Tq'"{t))v+ j Vg^(f) ■ Vu + [ X{t)q'^{t)v 
Jfi Jq Jfi 

- j dtqfit) vr + j Vrgr(t) • Vrfr + j Ar(t) qf{t) vr 
= / ipQ{t)v + / (^s(^)Pr for a.a. t G (0,T) and every (u,ur) ^ 


(fNq'^ + rg'^) (T) v+ QriT) ur = 0 for every {v, ur) £ • 


(2.40) 


(2.41) 

(2.42) 


More precisely, in [10] the above theorem is proved with the particular choice 

A = /"(j/^), Ar = /r(l/f), ifq = bQ{y^ - zq) and - ^s) (2.43) 

where y'^ and are the components of the state associated to an optimal control uf. 
However, the same proof is valid under assumption fl2.38p . 

Finally, imi Thm. 2.6] gives a necessary condition for uf to be an optimal control in 
terms of the solution to the above adjoint system corresponding to fl2.43p . This condition 
reads 

/ {Qt + boTif){vr — uf) > 0 for every vr G Had- (2.44) 

Jt. 

In this paper, we hrst show the existence of an optimal control ur- Namely, we prove 
the following result. 

Theorem 2.5. Assume fl2.3l) - fl2.7l) and fl2.11l) - fl2.13p . and let 3 and Had be defined by 
(I2.28P and (I2.29p under the assumptions (12.271) and (12.301) . Then there exists some ur G 
Had such that 

3{y, yr, ur) < 3{y, yr, Rr) for every ur G Had , (2.45) 

where y, yr: U ond yr are the components of the solutions {y,yr,w) and {y,yr,w) to the 
state system fl2.18p - 02.201) corresponding to the controls ur and ur, respectively. 
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Next, for every optimal control Mr, we derive a necessary optimality condition like 
(12.44 jl in terms of the solution of a generalized adjoint system. In order to make the last 
sentence precise, we introduce the spaces 

W ;= L2(0, T; Vn) C (if^(0, T; V*) x ll\0, T; Cr)) (2.46) 

Wo;={(M,Mr) e W: (m, Mr)(0) = (0, 0)} (2.47) 

and endow them with their natural topologies. Moreover, we denote by ((•, •)) the duality 
product between Wq and Wq. We will prove the following representation result for the 
elements of the dual space Wq. 

Proposition 2.6. A functional F : Wq —)■ M belongs to Wq if and only if there exist A 
and Ar satisfying 

Ae {H\0,T;V*)nL‘^{0,T;V)y and Ar G (i7^(0, T; P;) n ^^(0, T; W))* (2.48) 
{{F, {v, Mr))) = (A, v)q + (Ar, Mr)s for every {v, vr) e Wq , (2.49) 

where the duality products (•, ■)q and (•, • )e are related to the spaces X* and X with 
X = H\0,T-,V*)nL'^{0,T;V) and X = H\0,T-,V^) H L‘^{0,T-,Vr), respectively. 


However, this representation is not unique, since different pairs (A, Ar) satisfying fl2.48p 
could generate the same functional F through formula fl2.49p . 

At this point, we are ready to present our result on the necessary optimality conditions 
for the control problem related to the pure Cahn-Hilliard equations, i.e., the analogue of 
(I2.44P in terms of a solution to a generalized adjoint system. 

Theorem 2 . 7 . Assume (12.31) - (12.711 and (I2.11l) - (l2.13p . and let d and Had be defined by 
(I2.28P and (12.291) under the assumptions (I2.27P and (I2.30p . Moreover, let ur be any 
optimal control as in the statement of Theorem \2.5[ Then, there exist A and Ar satisfying 
(I2.48p . and a pair {q,qr) satisfying 


as well as 


q e L°°(0,T; W) C ^^(0, T; G) 

gr eL°°(0,T;i7r)nL2(0,T;Vr) 

(g, gr)(t) G Vq for every t G [0, T], 


{dtv{t),Xq{t)) dt + / {dtvr{t),qr{t))rdt 


(2.50) 

(2.51) 

(2.52) 


such that 


+ / Vg-VM+ / Vrgr • Vrfr + (A,m)q + (Ar,Mr); 


'Q 


= / / (pT.vr for every {v,vr) eWo, 


'Q 


(gr + boUr){vr - Mr) > 0 for every vr G Had- 


(2.53) 


(2.54) 


In particular, if bo > 0, then the optimal control ur is the -projection of —q^/bo onto 
Had- 
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One recognizes in fl2.53p a problem that is analogous to fl2.4ip - fl2.42p . Indeed, if A, Ar 
and the solution [q, q-p) were regular functions, then its strong form should contain both 
a generalized backward parabolic equation like fl2.4ip and a hnal condition for (!Ng, qp) of 
type fl2.42p . since the dehnition of Wq allows its elements to be free at t = T. However, 
the terms and Afgf are just replaced by the functionals A and Ar and cannot be 
identihed as products, unfortunately. 


3 Proofs 

In the whole section, we assume that all of the conditions fl2.3l) - fl2.7p and fl2.lip - fl2.12p 
on the structure and the initial datum of the state system, as well as assumptions fl2.27p 
and fl2.30l) that regard the cost functional fl2.28p and the control box fl2.29p . are satished. 
We start with an expected result. 

Proposition 3.1. Assume Up G Hp) and let {y'^,yp,w'^) be the solution to the 

problem fl2.18l) . fl2.20l) and fl2.23p associated to Up. If Up converges to up weakly in 


77^(0,T;77r) as r 

N 0, then 


y'" 

weakly star m 77^(0, T; V*) n L“(0, T; V) n T; 



and strongly in C°([0, T]; 77) fl 71^(0, T; V) 

(3.1) 

yf yr 

weakly star %n 77^(0, T; 77r) n 7:~(0, T; Hr) n 7:^(0, T; 772(r)) 



and strongly in C'°([0, T]; 77r) H L'^{0,T;Vp) 

(3.2) 

w'^ ^ w 

weakly star in L'^{0,T;V), 

(3.3) 


where {y,yp,w) is the solution to problem fl2.18l) - fl2.20p associated with up. 

Proof. The family {wf} is bounded in 77^(0, T; 77r). Thus, the solution {y'^,yf,w'^) satis- 
hes fl2.25p for some constant Cq, so that the weak or weak star convergence specihed in 
fl3.1D - fl3.3p holds for a subsequence. In particular, the Cauchy condition fl2.20p for y is 
satished. Moreover, we also have rdty'^ —?■ 0 strongly in L^(0,T;77) as well as t ^ 

and fp{yf) fp weakly in 77^(0, T ; 77) and in L^(0, T; Hp), respectively, for some f and .^r- 
Furthermore, y'^ and yf converge to their limits strongly in L^(0,T; 77) and 71^(0, T; 77r), 
respectively, thanks to the Aubin-Lions lemma (see, e.g., m Thm. 5.1, p. 58], which also 
implies a much better strong convergence [2T], Sect. 8, Cor. 4]). Now, as said in Section [21 
we can split /' as /' = /9 + tt, where (3 is monotone and vr is Lipschitz continuous. It 
follows that n{y'^) converges to 7i{y) strongly in L‘^{0,T; H), whence we obtain that also 
(3{y'^) converges to — 7i{y) weakly in L'^{0,T; H). Then, we infer that f — 7i{y) = (3{y) 
a.e. in Q, i.e., f = f'{y) a.e. in Q, with the help of standard monotonicity arguments 
(see, e.g., [H Lemma 1.3, p. 42]). Similarly, we have .^r = fp{yr)- Therefore, by starting 
from fl2.2ip and fl2.24p written with uf in place of up, we can pass to the limit and obtain 
fl2.2ip - fl2.22p associated to the limit control Mr- As the solution to the limit problem is 
unique, the whole family {y'^,yf,w'^) converges to {y,yp,w) in the sense of the statement 
and the proof is complete. □ 


Corollary 3.2. Estimate fl2.25p . written formally with r = 0, holds for the solution to 
the pure Cahn-Hilliard system fl2.18l) - fl2.20p . 
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Proof. By applying the above proposition with Wp = uy and using fl2.25p for the solution 
to the viscous problem, we immediately conclude the claim. □ 

Proof of Theorem 12.51 We use the direct method and start from a minimizing se¬ 
quence {wr.n}- Then, Mr,n remains bounded in H), whence we have Mr,n —t Mr 

weakly in Hy) for a subsequence. By Corollary 13.21 the sequence of the corre¬ 

sponding states {yn,yr,n,Wn) satishes the analogue of fl2.25p . Hence, by arguing as in the 
proof of Proposition 13.11 we infer that the solutions (i/„, |/r,n, Wn) converge in the proper 
topology to the solution {y,yY,w) associated to Mr- In particular, there holds the strong 
convergence specihed by the analogues of fl3.ip and fl3.2l) . Thus, by also owing to the 
semi continuity of 3 and the optimality of Mr,n; we have 

2{y,yY,UY) < hminfa(2/n,2/r,n,Rr.n) < 3(2/,2/r,Rr) 

n—>-OD 

for every Mr G Uad, where y and i/r are the components of the solution to the Cahn-Hilliard 
system associated with Mr. This means that Mr is an optimal control. □ 

Proof of Proposition 12.61 Assume that A and Ar satisfy fl2.48p . Then, formula fl2.49p 
actually defines a functional F on Wq. Clearly, F is linear. Moreover, we have, for every 
(M,Mr) G Wo, 


|(A,m)q -h (Ar,Mr)s| 

^ ll^ll(Rl(0,T;y*)nL2(0,T;V))* 11^||Ml(0,T;V*)nL2(o,T;y) 

+ ll^r||(Mi(o,r;yp*)nL2(o,T;yr))* ll'^l|Mi(o,r;yp*)nL2(o,r;yr) 

^ (l|l^|l(Rl(0,T;y*)nL2(0,T;y))* + ||Ar||(Ml(0,T;yp*)nL2(0,T;yr))*) ll('^Wr)||w , 

so that F is continuous. Conversely, assume that F G Wq. As Wq is a (closed) subspace 
of W := T- V*) n L2 (o, T; V)) x T; Vp*) C ^^(o, T; Vr)), we can extend F to 

a linear continuous functional F on W. Then, there exist A and Ar (take A(m) := F{v, 0) 
and Ar(Mr) := .F(0,Mr)) satisfying fl2.48p such that 

{F, (m, Mr)) = (A, v)q + (Ar, Mr)s for every (m, Mr) G W , 

where the duality product on the left-hand side refers to the spaces (W)* and W. Since 
{{F, (M,Mr))) = {F, (M,Mr)) for every (M,Mr) G Wq, fl2.49D immediately follows. □ 

The rest of this section is devoted to the proof of Theorem 12.71 on the necessary 
optimality conditions. Therefore, besides the general assumptions, we also suppose that 

Mr is any optimal control as in Theorem 12.51 (3-4) 

that is, an arbitrary optimal control Mr is fixed once and for all. In order to arrive at the 
desired necessary optimality condition for Mr, we follow [2] and introduce the modified 
functional 3 dehned by 

3iy,yr,UY) := 3iy,yr,UY) + - ||Mr - Rr|li2(s) • (3.5) 

Then the analogue of Theorem 12.31 holds, and we have: 
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Boundary control problem for the pure Cahn-Hilliard equation 


Theorem 3.3. For every r > 0, there exists some Wp G Had such that 

Siv'", vh “r) < Vr, Ur) for every ur G Had , (3.6) 

where y'^, y^, y'^ and yf are the components of the solutions and {y'^,yf,w'^) 

to the state system (I2.18p . (12.201) and (I2.23p corresponding to the controls wf and ur, 
respectively. 


For the reader’s convenience, we fix the notation just used and introduce a new one 
(which was already used with a different meaning earlier in this paper): 

is an optimal control as in Theorem 13.31 (3.7) 

{y'^,y^,w'^) is the solution to 02.181) . O2.20p and 02.231) corresponding to Up (3.8) 
(y'^,yf,TiJ'^) is the solution to 02.18p . O2.20p and 02.231) corresponding tour- (3.9) 

The next step consists in writing the proper adjoint system and the corresponding neces¬ 
sary optimality condition, which can be done by repeating the argument of HD]. However, 
instead of just stating the corresponding result, we spend some words that can help the 
reader. The optimality variational inequality is derived as a condition on the Frechet 
derivative of the map (defined in a proper functional framework) up eA 3(y, yr, up), where 
the pair (y, yr) is subjected to the state system. Thus, this derivative depends on the 
Frechet derivative of the functional (y,yr,ur) eA 3(y,ur,ur), which is given by 

[D3(y,yr,ur)](k,kr,hr)] = bQ / (y-ZQ)k + bj: / (yr-ZE)kr+ / (boUp + (up - ur)} hr . 

Jq Jy. Jy 


Hence, the argument for 3 differs from the one for 3 only in relation to the last integral. 
In other words, we just have to replace bQUp by bQUp + {up — Tip) in the whole argument 
of [To]. In particular, the adjoint system remains unchanged. Here is the conclusion. 

Proposition 3.4. With the notations (13.71) - (13.81) . we have 

[ {qf + bouf -f (wp - Fr)) (vp - uf) > 0 for every vp G Had , (3.10) 

Jy 

where gf is the component of the solution (g^,gf) (I2.40p - (l2.42p corresponding to up = 

Up with the choices X = A"", Ar = Ap, (pq = pq and = Ts specified by 

y = rm, Af = = 6g(r - ^q) “^i = 65,(9? - ze). (s.h) 


Thus, our project for the proof of Theorem 12.71 is the following: we take the limit in 
fl3.10p and in the adjoint system mentioned in the previous statement as r tends to zero. 
This will lead to the desired necessary optimality condition fl2.54p provided that we prove 
that the optimal controls uf converge to Tip. The details of this project are the following. 

i) There hold 


Up Up 

¥ 


(3.12) 


weakly star in Hp) nT°°(S) and strongly in T^(S) 

weakly star in 77^(0, T; V*) n T°°(0, T; V) n T; 772(H)) 
and strongly in C°([0, T]; 77) D T2(0, T; V) 


(3.13) 
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yl weakly star in T; Hr) n L“(0, T; Vr) n ^^(0, T; H^{T)) 

and strongly in C°([0, T]; Hr) C T; Vr) (3.14) 

w'^ weakly star in L^(0,T; V) (3.15) 

^ q weakly star in L°°(0,T; C*) n L^(0,T; G) (3.16) 

^ Q'r weakly star in L°°(0, T; ilTr) n L^(0, T; Cr), (3-17) 

as well as 

H(y^, y?, ul) diV, Vr, Ur), (3.18) 


at least for a subsequence, and {y,yr,w) solves problem (I2.18p - (l2.20j) with Mr = Rr- 

ii) The functionals associated with the pair (A'^g'^, Af, gf) through Proposition [2]6] converge 
to some functional weakly in Wq, at least for a subsequence, and we then represent the 
limit by some pair (A, Ar), so that we have 

(A^g^, v)q + (Apgr, Mr)s ^ (A, v)q + (Ap, Mr)E for every (n, vr) e Wq. (3.19) 

Hi) With such a choice of (A, Ap), the pair (g, gr) solves fl2.52l) - fl2.53p . 
iv) Condition fl2.54l) holds. 

The main tool is proving a priori estimates. To this concern, we use the following rule 
to denote constants in order to avoid a boring notation. The small-case symbol c stands 
for different constants that neither depend on r nor on the functions whose norm we want 
to estimate. Hence, the meaning of c might change from line to line and even in the same 
chain of equalities or inequalities. Similarly, a symbol like cs denotes different constants 
that depend on the parameter 6, in addition. 

First a priori estimate. As wf G IXad and Theorem 12.21 holds, we have 

\\Ur\\m{0,T-,Hr) + ll^’'l|Rl(0,r;V*)nI,°°(0,T;V)nL2(0,r;R2(0)) 

+ ll^rllRh0,T;Rr)nL°°(0,T;Vr)nL2(0,T;R2(r)) + \\u!^\\l^{0,T-,V) 

+ II/^(y’') 11^2(0,T;R) + \\fr{yr)\\L^(0,T-,Hr) + '^^'^^ |||| L2 (o,T;R) < C. (3.20) 


Second a priori estimate. For the reader’s convenience, we explicitly write the adjoint 
system mentioned in Proposition 13.41 as well as the regularity of its solution, 

q^ e H\0,T;H)nL'^{0,T;H\n)), q^. e H\0,T; Hr) f] L'^{0,T; H\T)) (3.21) 

(?^,5r)(s) e Vn for every s e [0,T] (3.22) 


- / ai(K(g"(s)) + rg"(s))M+ / Vg"(s) • Vm + / A"(s)g"(s)M 

Ju Jq Jn 

- j dtq^s) ur + J Vrg?(s) • VrMr + j Xr{s)qr{s)vr 

= / (Pq{s)v + / <yC 2 (^)w for a.e. s G (0,T) and every (M,Mr) £ A’o (3.23) 

Jn Jr 

where A^ =/"(^^), Af =/f'(^f), = bqiy^ - zq) and = bj:{y^ - zj:) 

/ (Afg’" -I- rg’") (T) v+ qriT) Mr = 0 for every (m, vr) G Vf 


'n- 


(3.24) 
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Now, we choose v = g^(s) and Pr = Q'f('S)) and integrate over (t,T) with respect to s. 
Recalling fl2.35p and now reading := Q x {t,T) and E* := F x {t,T), we have 


+ l [ k'WlY [ |v«T+ f 

^ ^ Jn JQt JQt 

+ A / kf(^)r+ f |Vrgfr+ f -^rkfl^ 

" «/r St St 





<1^+1 

Jt. JT.t 


< 




IgpP + c 


(3.25) 


where the last inequality follows from fl3.20l) . By accounting for fl2.5l) . we also have 


>-c \q 


r|2 


and 


Arlo'rl ^ c / |gp 


r |2 


^ Qt 


^ Qt 


' Qt 


'St 


We treat the volume integral (and the same on the right-hand side of fl3.25lF invoking the 
compact embedding V C H. We have 




|Vn|^ cs\\v 


for every n G R and 5 > 0. 


Hence, we deduce that 


/ \Vq-\^ + cs rU{s)\\lds. 

Qt J Qt J t 

Therefore, by combining, choosing 6 small enough and applying the backward Gronwall 
lemma, we conclude that 


lk’"l|L°°(0,T;U*)nL2(0,r;U) + lkrl|L°°(0,T;Rr)nL2(0,T;Ur) +< C. (3.26) 


Third a priori estimate. Take an arbitrary pair (n, up) G iF^(0, T; TCq) nL^(0, T; Vf^), 
and test fl3.23p by v{s) and nr(s). Then, we sum over s G (0,T) and integrate by parts 
with the help of 03.241) . so that no integral related to the time T appears. In particular, 
if (n,nr) £ AVq, even the terms evaluated at t = 0 vanish and we obtain that 


(Ng'^ -I- Tq'^)dtV -|- / Vg'^ ■ Vn -|- / X^q'^v + / gp^iUp -|- / Vgp • Vfp -|- / A'^gpUp 


'Q 


'Q 


>Q 


= / / <Fs^r. 


(3.27) 


>Q 


Therefore, we have, for every (n,np) G Wq, 


[ A"g"n+ [ A'gfnp 
Jq Jt, 

< ||Ng'^ -I- rg'^||p2(o,T;y) ||^tP|U2(o,T;y*) + |k’'||L2(o,T;V) I1 p|U2(o,T;U) 
+ Ikf llL2(0,T;yr) ll^tW||L2(o,r;Vp*) + \\qr\\L^(0,T-,Vr) llw||L2(o,r;Vr) 
+ 117^011^2(0^;^) II^I|l 2(0,T;R) + 117^211^2(0,T;Rr) l|w||L2(o,T;Rr) ' 
















Colli 


Gilardi 


Sprekels 


15 


Now, by assuming r < 1, we have ||!Nu + tv\\v < c||u||* + t||p|| v < c||u||\/ for every v eV 
with zero mean value (see fl2.36p ). Therefore, by accounting for fl3.20p and fl3.26p . we 
conclude that 


/ yq^v+ / yqlw 

IQ JS 


< c||(u,ur)||w for every (u,ur) e Wq. 


(3.28) 


Conclusion of the proof of Theorem 12.71 From the above estimates, we infer that 


wf Ur 

weakly star in T; Hr) 0 L°°(S) 

(3.29) 

¥ 

weakly star in 77^(0, T; V*) 0 L°°(0, T; V) 0 ¥{0, T; H^{n)) 



and strongly in C°([0, T]; H) fl L^(0, T; V) 

(3.30) 

¥ hr 

weakly star in 77^(0, T; 77r) n L“(0, T; Vr) n T; H^{T)) 



and strongly in C°([0, T]; 77r) nL^(0,T;l/r) 

(3.31) 

w'^ ^ w 

weakly star in L^(0, T; V) 

(3.32) 

q'^ ^ q 

weakly star in L°°(0, T; V*) n T; V) 

(3.33) 

Qr Qr 

weakly star in L°°{0,T; Hr) D L‘^{0,T;Vr) 

(3.34) 

rq'^ —)■ 0 

strongly in L°°(0, T; H) 

(3.35) 


at least for a subsequence, and {y, yr, w) is the solution to the problem fl2.14p - fl2.20p 
corresponding to ur, thanks to Proposition 13.11 Notice that fl3.33p - fl3.34p coincide with 
fl3.16p - fl3.17p and that fl3.12p - fl3.15p hold once we prove that ur = ur and that uf con¬ 
verges strongly in L^(S). 

To this end, we recall the notations fl3.7p - fl3.9l) . and it is understood that all the limits 
we write are referred to the selected subsequence. By optimality, we have 

3{y,yr,ur) < d{y,yr,ur) and < a(y^,yf,Mr)- 

On the other hand, fl3.29p - fl3.3ip and Proposition 13.11 applied with uf = ur yield 

3(|/,2/r,Rr) < liminf3(^^,^r,Mr) and \im3iy'",yr,ur) = 8{y,yr,ur). 

By combining, we deduce that 

8{y,yryUr) + - llwr - Wr||i 2 (s) < 8{y,yr,ur) + - Hur - wr||L 2 (s) 

= 3{y,yr,ur) < hminfa(y^,y?,Mr) < limsup hr) 

< hmsupa(y^,yf,Mr) = limsup a(y^, yf, hr) = 3{y,yr^ur). 

By comparing the hrst and last terms of this chain, we infer that the L^(S)-norm of 
Mr — Ur vanishes, whence ur = ur, as desired. In order to prove the strong convergence 
mentioned in fl3.12p . we observe that the above argument also shows that 

hminf3(^'^,yf,hp) = limsup a(y^, hr) = 3iy,yr,ur). 

Notice that this coincides with fl3.18p . From the strong convergence given by fl3.13p and 
fl3.14p . and by comparison, we deduce that 
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whence also 


lim sup ^ 


\u 


t|2 


< lim sup ( — 


\u 


t|2 




\Ut 




bo 

2 


Rp S 


^ < lim inf 


b_o 

2 


\uy\^. 


Therefore, we have 




whence 



0 , 


and (I3.12I) - (I3.15I) are completely proved. 

Now, we deal with the limit {q,qr) given by fl3.33p - fl3.34p . i.e., fl3.16p - fl3.17p . Clearly, 
fl2.52p holds as well. Furthermore, as ||!Np*||y < c||t*||* for every n* G V* with zero mean 
value (see 02.3311 ). and since the convergence fl3.35p holds, we also have 


—)■ Ng weakly star in L°°{0,T; H). 


On the other hand, fl3.28p implies that the functionals F’" G Wq dehned by 

{{F% (p, nr))) := (A^g", r)q + (A^gf, Rr)^ , 

i.e., the functionals associated with (A'^g'^, Af, gf) as in Proposition l2.61 are bounded in Wq. 
Therefore, for a subsequence, we have —)■ F weakly star in Wg, where F is some element 
of Wq. Hence, if we represent F as stated in Proposition 12.61 we hnd A and Ap satisfying 
fl2.48p and fl3.19p . At this point, it is straightforward to pass to the limit in fl3.27ll and in 
fl3.1Up to obtain both fl2.53p and fl2.54p . This completes the proof of Theorem 12.71 □ 

Remark 3.5. The above proof can be repeated without any change starting from any 
sequence n 0. By doing that, we obtain: there exists a subsequence such 

that fl3.12p - fl3.18p hold along the selected subsequence. As the limits ur, y, y^, w and 
3{y, Vr^ Rp) are always the same, this proves that in fact fl3.1211 - 03.1511 as well as fl3.18p hold 
for the whole family. On the contrary, the limits g and gp might depend on the selected 
subsequence since no uniqueness result for the adjoint problem is known. Nevertheless, 
the necessary optimality condition fl2.54p holds for every solution (g, gp) to the adjoint 
problem that can be found as a limit of pairs (g"^, gf ) as specihed in the above proof. 
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